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We study screening properties of the two-dimensional electron gas with Rashba spin-orbit cou- 
pling. Calculating the dielectric function within the random phase approximation, we describe the 
new features of screening induced by spin-orbit coupling, which are the extension of the region of 
particle-hole excitations and the spin-orbit-induced suppression of collective modes. The required 
polarization operator is calculated in an analytic form without any approximations. Carefully de- 
riving its static limit, we prove the absence of a small-g anomaly at zero frequency. On the basis 
of our results at finite frequencies we establish the new boundaries of the particle-hole continuum 
and calculate the SO-induced lifetime of collective modes such as plasmons and longitudinal optical 
phonons. According to our estimates, these effects can be resolved in inelastic Raman scattering. 
We evaluate the experimentally measurable dynamic structure factor and establish the range of 
parameters where the described phenomena are mostly pronounced. 
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I. INTRODUCTION 

Fundamental issues of interaction effects in the two- 
dimensional electron gas (2DEG) are at the center of 
discussions since the early days of its fabrication^. Quite 
generally, screening described by the dielectric function 
forms the basis for understanding a variety of static 
and dynamic many-body effects in electron systemsS. 
The dielectric function of the 2DEG was computed a 
long time ago by Stern^ within the random phase ap- 
proximation (RPA). Different quasiparticle and collec- 
tive (plasma) properties deduced from those expressions 
were confirmed experimentally soon after—. Recent ex- 
periments measuring plasmon dispersion, retardation ef- 
fects, and dampin g-'-'^ unambiguously show the impor- 
tance of correlations between electrons. 

More recently, the possibility of manipulating spin in 
2DEG by nonmagnetic means has generated a lot of 
activity'*. The key ingredient is the Rashba spin-orbit 
(SO) coupling^ tunable by an applied electric fieldiS. A 
recent example of SO-induced phenomena which has at- 
tracted much attention is the spin-Hall effect^!. 

In this context the study of interplay between electron- 
electron correlations and SO coupling in 2DEG be- 
comes an important problem. In the preceding pa- 
pers it has been already discussed how SO coupling 
affects static screeningiS, plasmon dispersion and its 
attenuationiii^ii^, and Fermi-liquid parametersi^. How- 
ever, the approaches of these papers as well as of many 
other papers on Rashba spin-orbit coupling are based on 
various approximations. We can outline the most pop- 
ular ones: a linearization of spectrum, also known as ^- 
approximation; an expansion in SO coupling parameter 
up to the lowest nonvanishing contribution; reshuffling 
the order of momentum integration and evaluation of the 
zero-frequency limit in the polarization operator; a com- 
bination of any of those approximations. In our opinion, 
their accuracy is not comprehensively discussed. Since 



its lack might seriously affect a description of SO-induced 
phenomena, this issue should be thoroughly investigated. 

In the present paper, we study the effects of the dy- 
namic screening in 2DEG with Rashba SO coupling de- 
scribed by the RPA dielectric function in the whole range 
of momenta and frequencies. The main part of our pa- 
per is devoted to the analytic evaluation of the polar- 
ization operator, which does not employ any approxima- 
tions. The knowledge of the dielectric function allows 
us to predict new features in the directly observable dy- 
namic structure factor. In particular, we observe a SO- 
induced extension of the particle-hole excitation region 
and calculate the SO-induced broadening of the collec- 
tive excitations such as plasmons and longitudinal opti- 
cal (LO) phonons. We obtain the values of lifetimes that 
lie in the range of parameters experimentally accessible 
by now in the inelastic Raman scattering measurements. 

On the basis of our analytic results, we also revisit 
the earlier approaches, estimating their accuracy and es- 
tablishing the limits of their applicability. Although the 
approximations mentioned above usually work well for a 
conventional 2DEG (without SO coupling), we demon- 
strate that in case of the Rashba spectrum they should 
be applied with caution. In our paper we discuss the 
subtle features of the 2DEG with a SO coupling warning 
about this. Special attention is focused on a derivation of 
a zero-frequency (static) limit of the polarization opera- 
tor. We present an analytic result, which, however, does 
not contain an anomaly at small momenta predicted in 
Ref.lli 

The paper is organized as follows. In Sec.^we intro- 
duce the main definitions and notations. In Sec. IIIII we 
present a detailed analytic calculation of the polarization 
operator I{{q,Lo) and discuss the most important modi- 
fications generated by SO coupling. We also derive an 
asymptotic value of ImH(g,ti;) at small q and compare 
it to the approximate expressions known befor o^^i^^ . In 
Sec. II VI we implement a thorough analytic derivation of 
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the static limit \im^^()Il{q,u}) and make a conclusion 
about the absence of an anomaly at small momenta. In 
Sec. we study the directly observable dynamic struc- 
ture factor and evaluate SO-induced plasmon broaden- 
ing. The latter effect also causes a modification of the 
energy-loss function that is estimated as well. In Sec. IVII 
we calculate the lifetime of LO phonons generated by SO 
coupling. 



II. BASIC DEFINITIONS 

We consider a 2DEG with SO coupling of the Rashba 
typeS described by the single-particle Hamiltonian, 



III. POLARIZATION OPERATOR AT 
ARBITRARY FREQUENCY 

A. Evaluation of the polarization operator 

For a calculation of (jSJ it is useful to shift k ^ k — q in 
those terms of (O that contain np(e^'^^). At the same 
time we can shift the integration angle — > -I- tt in 
the same terms due to the momentum isotropy of the 
spectrum Conveniently regrouping all contributions 
in Eq. Q), we cast it into the form 



n„ 



A' 



(8) 



H 



2m* 



aii}\{a X k), 



(1) 



where n is a unit vector normal to the plane of 2DEG 
and h ~ 1. The dispersion relation is SO split into two 
subbands labeled by /i = ±, 

-+HaRk= — (2) 



2m 



2m* 



2m* 



These subbands have the distinct Fermi momenta fc^ = 
hp — fikji and the same Fermi velocity vp — kp/m*. 



Here we denote kn = m*aji and kp = ^2m*Ep 4- 
and assume that 2kpi < kp. 

The effective Coulomb interaction is V^^^ — Vqle^^, 
where Vq = 27re^/((7£oo), and £oo is the (high-frequency) 
dielectric constant of medium. The dielectric func- 
tion Squj = Ree^„ -I- jlmSgi^ describes effects of dy- 
namic screening, and in the random phase approximation 
(RPA) it is given by^ 

(3) 



where n^j^ is a polarization operator. In the presence of 
SO coupling, the latter is a sum 



E 

X,M=± 



/ J ^^quj 



of contributions 



m.:^ = lim 



d^k np[^^)-np{e^^^) 



5^0 7 (27r)2 t^ + iJ + eJ^ 



k,k+q 



(4) 



, (5) 



where the indices x = + and x = ^ correspond to the 
intersubband and intrasubband transitions, respectively. 
The form factors 

•^k.k+q = ^ [1 + X C0s((/)k - (/)k+q)] (6) 

originate from the rotation to the eigenvector basis, and 
^ . Ik|+x|q| 



COS((/)k 



!^k+q^ 



|k + q| 

X — COs(0k ~ 0q) = COS ( 



(7) 



Throughout the paper we will use the dimensionless 
units y = kn/kp, z = q/2kp, v = k/kp, and w = 
m*LLj/2kp. 



where 

n^:, , , = lim 



(9) 



•^k,k+q 



k,k+q 



'"k+q 



+ X{UJ + iS) 



*^k fk+q 



+ X{uj + iS) 



and the index A effectively labels the contributions from 
the different ("in" and "out") Fermi functions. 
In the limit of zero temperature we obtain 

1 + COs((?ik - </>k-|-q) , 1 - COs((/)k - 0k+q) 



'"k+q 



-I- A(w + iS) 



"^k "^k+q 



+ X{uj + iS) 



After the intermediate steps (jAl|) , (jA2|) we cast (|1U|I into 
the form 



-ImH'f^ 



qu^X 



vg^{v,z,w,y)dv, (11) 







-ReU 



'jMv,z,w,y)dv, (12) 



where v = i/2d = m* /{2t:) is the density of states in 
2DEG per each spin component; the indices i = 1,2,3,4 
correspond to {/i, A} = {-,+}, +},{-,-}, -}, 
and the functions gi, fi are defined by 



sign(2wza; — fiyv + 2(z^ — Aw)) 



X (x + Si) S{x'^ + PiX + 7i), 

a , x + s,. 



27r Jq x^ + f3iX + ji ' 



(13) 
(14) 



with the coefficients 
V- fiy 



a = 

7» = 



(z^ — Xw) — ^yv 



2v'^z ' * {v — IJ'y)z 

2(z2 - Xw) - fiy{v + fiy) 



{z^ — Xw)'^ — ^yv{z^ — Xw) — z'^y'^ 



(15) 
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In (|13|l the Dirac delta function is denoted by (5 (...), as 
usual. 

The equation + fiiX + 7i = has the roots labeled 
by cr = ±, 



2wz 



(16) 

The roots Xia and A2(t are always real, while and A4ct 
become complex in the ranges y — 2y/w < v < y + 
and —y — 2y/w < v < — y + 2y/w, respectively. 

All integrals in and (|12|l receive a contribution 
from those ranges of v where Xia are real. Additionally, 
the integrals J vf^^^dv receive a contribution from the 
ranges with complex Xs^ia-- In order to take into account 
both types of contribution, we represent the function fi = 
fi + fi' as a sum of // = fMPf - and //^ = 
/.e(47, - (3f). Obviously, ///a = 0. 

Let us also introduce Rell^ and RcII^^ that are ob- 
tained by integrating the functions J2i "^fi and ^ ■ w//^, 
respectively. Thus, the full real part of the polarization 
operator is given by the sum of the two, i.e.. 



Ren = Ren' + Ren 



(17) 



Note that in the limit w — > all the roots Xia become 
real, and it might seem that Ren^^ vanishes in the static 
limit. However, it is not the case, and Ren-^^ does give a 
finite contribution as w — > 0. We will discuss this point 
in a deep detail in the next section. 
For (3f > 47i we have 



c/ 9 ^ N S(x — Xi+) + S(x — Xi^) , 
(5(x2 + fi,x + 7,) = ^ p ^,(18) 



1 



+ /3iX + 7j Ai+ - Ai_ \x- A, 



Ai+ — Ai_ 
Aj+ + 6i Aj_ + 5j 



X - Aj_ 



(19) 



Using the integrals (jA8|) and \K%\ we establish that 
Hadia, fi = E^/ii' and 



XC, 



9i<j = 



Xia + Si 



:e(i- |A„ 



(20) 



X Q{{v + ^.yf + 4Aw) sign(?; + cr^/{v + ^y)^ + AXw), 



fi = 



C, 



Xia + Si 



Ai+ - Xi- 
For Pf < 4'-fi we have 



1 



e(|A,,| - l)sign(A,,). (21) 



1_ 



--Ren" = 



vf^'dv+ I vfi'dv, (22) 



where I?zp are defined in the following way: 



y > 4w 
Aw > y"^ n 1 > Aw 
Aw > 1 



V_ = [y-2V^,y + 2y^,{23) 
V^^[0,y + 2V^l (24) 
l?_ = [0,l + y], (25) 



and 

y2 > 4w : V+ = II), (26) 
V+ = [0,-y + 2V^l (27) 
v+ = [o,i-y]. (28) 
Applying the table integral (|A10|) . we obtain 

1 



Aw > y ni > Aw 
Aw > 1 



f^i{v,z,w,y) 



2z^{v^^~v){v~vi^) 



(29) 



where 



= —^y ± (z + w/z), 



"3,4 = 



±z 



tiy±z' 



(30) 

(31) 
(32) 

(33) 



and the indices /i = — and /i = + correspond to /g and 
fl', respectively. 

The numerical evaluation of l|22|l can be performed 
with a controlled and sufficiently high accuracy, since the 
functions (|29() may diverge only near the actual integra- 
tion edges —iJ,y±2y/w. On the other hand, one can try to 
find an analytic expression for (|22|l . An alternative rep- 
resentation of /1 4 , which is equivalent to (|29|l and more 
suitable for further analytic evaluation, is introduced in 
Appendix ^ 

As for Imn = ^ i^X]i=i/o^ vgidv and RcH^ = 
— i/^^^j^ '^^vfidv, their analytic evaluation is pre- 
sented below. Making a change of variables, 



1 



r = — 
2 



-pL{v + fiy) + a^/(v+Jly)^ + AXw 
we establish the relations ljA3|l - (jA7p . and thus deduce 



(34) 



1 

-Imn = VfT / dTC+{T) (35) 

- 9(1 -4u;)^ / dT£-{T) 
f'<-+-(y) 

+ 2Q{y'^~Aw) / dTC-{T), 



and 



Ren^ = J2 



+ {y) 

Qil-Aw)J2 



dTTZ+ir) 



(36) 



a,fj, 



dTn-{T) 

2Q{y^-Aw)y] [ drn-ir), 



Air++(0) 
.T-+(0) 
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FIG. 1: The domains A, B, C, V [see Eqs.lgSJ-lllEU corre- 
sponding to the different orderings of the roots I37I I. 



where 



n,2 

£(r) = 
7^(r) = 



±w/z, T3^4 

— —X + <J\/ + AXw 

C{t) sign(T'^ + ± w), 
'R{t) sign(r2; =F w{t + y)/ z) 

1 (t-t-3)(t-T4) 

"'I 



1 (t^T3)(t-T4) 




k=l 



In Fig. n the quarter-plane (z > 0, w > 0) is di- 
vided into the domains A — {{z,w)\w < z{z — y)}, 
B = {{z,w)\'w > z{z — y)nw < z{z + y) Dw > z{y — z)}, 
C = {{z,w)\w > z{z + y)}, V = {{z,w)\w < z{y - z)}, 
which are specified by an ordering of the roots Tk (|37() : 



A 
B 
C 
V 



T4 < T2 < Tl < T3, 

T4 < T2 < T3 < Tl, 

T2 < T4 < T3 < Tl, 

T4 < T3 < T2 < Tl. 



(43) 
(44) 
(45) 
(46) 



In each domain one should compare t^ and Ta-\{x) in 
order to establish actual limits of integration in (|35|l and 
(jSnj. After that, it becomes possible to write down Imll 
and Rell-^ in an explicit form. We refer to the Appendix 
where the necessary expressions for establishing the 
explicit form of (|35|l are presented. Similar expressions 
for H3t)|) can be easily derived from Ref. IT^ 



B. SO-induced extension of the particle-hole 
excitation region 

Let us analyze the most important modifications to the 
polarization operator induced by SO coupling. First of 



all, we are interested in establishing the new boundaries 
of a particle- hole continuum (or Landau damping region) , 
which is defined by the condition Imll ^ 0. They can 
be determined from a simple consideration of extremes 
of the denominators in Q. On the basis of this purely 
kinematic argument, it is easy to establish that due to 
SO coupling there appears a new wedge-shaped region 
of damping (shown in Fig. [21) . It is bounded by the two 
parabolas —(z — y)^ — (z — y) = W4{z) <w< wi{z) — {z+ 
y)^ + {z + y) and attached to the boundary wo{z) = z'^ + z 
of the conventional particle-hole continuum (obtained in 
the absence of SO coupling according to Ref.0). Another 



boundary w ~ z 



of the latter transforms into 



(z — j/)^ — [z — y) for nonzero y (this occurs at z > 1, and 
therefore it is not shown in Fig.EJ. 

An extension of the particle-hole continuum reflects an 
opened possibility for transitions between SO-split sub- 
bands. Therefore one can expect new SO-induced ef- 
fects of damping of various collective excitations (plas- 
mons, LO phonons) in the new regions of damping. 
These issues will be discussed in the subsequent sections. 
Below we present Imll in explicit form for the values 
w > wa{z), i.e., above the boundary woiz) of the conven- 
tional particle-hole continuum. Since u'o(z) > z^ -t- yz, 
we deal with the case of the roots' ordering H45|l corre- 
sponding to the domain C. Assuming that y does not 




0.02 



0.06 



0.1 z 



FIG. 2: SO-induced extension of the particle-hole continuum 
shown for y = 0.1. It is bounded from above by parabola 
'wi{z) [Eq. 14911 1 and bounded from below by the parabola 
W4{z) [Eq. II52II 1 and by the parabola wo{z) — z'^ + z oi the con- 
ventional particle-hole (Landau damping) region. The small 
darkened triangle indicates the region where the approxima- 
tion a{z,w) a{0,w) is applicable. 
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exceed the value ^ ^ ~ 0.3, we deduce from (|^ . 

1 

- -ImU = -e{w2iz) - 10)0(10 - Wi{z)) dr£(r) 

„i[l_7T+4in] 

- e(u;i(z) -'u;)e(w- W2(^)) / dT/:(T) (47) 

/.i[V^^-i] 

+ 6(^3(2) - w)e(w - 'u;4(z)) / dT/:(r). 

Applying the table integral ICSp for xx = Ti, Xa = X2 = 
T3, Xb — X3 = Ti, and X4 = T2, we obtain the analytic 
expression 

Imll = —Q{w2{z) — w)Q{w — W4{z))A{z — y) 



e(wi(z) - w)e{w - W2iz))A{-[l - VI + 4w])(48) 



+ e(u;3(^) - w)eiw - W4{z))A{-[VT^^ - 1]), 



where 



wi{z) = (z + yf + {z + y), (49) 

W2iz) = {z-yf-{z-y), (50) 

w:i{z) ^ -{z + yf + {z + y), (51) 

W4{z) = -{z-yf^{z-y), (52) 



and 



2z V w + zx 4zynv 
X [((?«/z-2;)2-z2)F(¥.(a;),fc) 

-((z«/z + z)2-y2)£;(^(^)^^) 

+2y(w;/z-y-z)n((p(x),n,fc)]. (53) 

The argument f{x) and the parameters k and n of the 
elliptic functions F, E and 11 are defined in the following 
way: 



, , / (a; + z + w)(w/z — w + z) 
^{x) = arcsinJ^ / , ^ (54) 



2(za; + w) ' 
2./u; 2z 



■\/(w/z + z)2 — t/2 



,(55) 



Note that (/?(— z — y) = and Lp{z — y) = -I, and hence 
A{z — y) is expressed in terms of the complete elliptic 
integrals F{§,k), E{^,k), and n(f ,n, fc). 



(O , and therefore leads to the kinematic extension of the 
conventional particle-hole continuum to a strip y — < 
w < y + y'^ parallel to the z axis. This approximation 
can be effectively expressed in the form 



U{q,uj) « -i-^a{0,uj), 



and considered as stemming from the identity 



n(<z,^) 



'«^o-('7,^), 



(56) 



(57) 



with aiji^uj) replaced by CT(0,tt'). The optical conductiv- 
ity <7{Q,uj) = <7^ can be easily foundi^. For example, its 
real part equals 



lOTT 



Thus, Eq. linnj yields 
1 



-imn 



(58) 



(59) 



We would like to argue that the approximation H59() 
works well only in a quite small region restricted by 
the conditions 103 (z) < w < W2{z) (see Fig. EJ. Since 
the parabolas 'W2{z) and u'3(z) intersect at z = i[l — 

\/l — 4?/2] « j/^, this region is represented by a small tri- 
angle between the points (0, y — j/^), (0, y -1- y^), (y^, y). 
Inside this triangle —IvaR/v is entirely determined by 
—A{z — y). Observing that fc^ « 2n « 4z^/w for small 
z < y^ and expanding the complete elliptic integrals with 
respect to z, we find an asymptotic value, 



i , , TTZ 



(60) 



which coincides with H59|l . except for the domain of ap- 
plicability. 

In Fig. O we compare the curves 167rRe[(T(zc, ?i')]/e^ 
calculated at the different values of Zc and y = 0.1, and 
plotted as a function of w. The unit-step between y — 
y^ = 0.09 and y + y'^ — 1.01 corresponds to the optical 
conductivity (|58|l . which is recovered at Zc = 0. One can 
observe that for z ^ y^ the approximation (t(z, w) ~ 
(t(0, w) leading to (|59|l works quite well, while for z^ y^ 
it completely breaks down. 



IV. STATIC LIMIT OF THE POLARIZATION 
OPERATOR 



C. Comparison with the approximate result 



Careful derivation of the static limit 



It is also worthwhile to compare our exact result 
with an approximation for small z commonly used in the 
literature (cf., e.g., Refs.0Q|l^. It neglects the square 
of the transferred momentum 5^ in the denominator of 



In this section we derive the static limit 

lim Ren(z,w) = lim Ren^(z,K7) + lim Ren^^(z,w). 

w — >0 w — >0 w — ^0 

(61) 
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167rRea(z=z^,w)/e^ 
1 

0.8 
0.6 
0.4 
0.2 



A 



III \ '.1 



0.05 



0.1 



0.15 w 



FIG. 3: The real part of conductivity for the values of 
Zc = 0.0, 0.025, 0.005, 0.01, 0.02, 0.03, 0.04 calculated at 
y = 0.1. The unit step corresponds to Zc — 0. The dashed 
line corresponds to Zc = t/^ = 0.01, which is the limiting value 
for the applicability of the approximation a{z,w) ~ a{0,w). 



After Eq. (|17|l we have already made an observation 
that Hni^_>o Ren-^^(z, w) might seem to vanish. Below 
we prove that, in fact, it does not vanish, but rather 
produces a finite and a very important contribution to 
lim^j^o Ren(z, w) for small z < y. 

Let us first identify the quantity lim^^o Ren^(z,'u;). 
One can see that if we put w = in all the roots 

Xia become real for all values of the integration variable 
V. Tracing back the derivation of the expression H14|) . it 
is easy to see that lim^j^o Ren^(z, w) is obtained, in fact, 
from the definition (cf. Ref. fl^ 



(2^) 



■^k,k+q- 

(62) 

We would like to emphasize that, in general, Tlq,uj+iS=o 
is not always the same as lim^^^o ^quj , and there might 
occur a specific situation when 



lini Hgj^ 7^ n<},u;+i(5=o- 



(63) 



Once lim„,_+o Ren^^(z, w) is nonzero, the subtle property 
(inSl, for example, holds for the 2DEG with Rashba SO 
coupling described by the Hamiltonian 

The correct static limit is, of course, given by 
limi^^ollgw- However, it is very tempting to put lo — 
from the very beginning. In principle, one can do that. 
But then, in order to protect oneself from a possible mis- 
take, it is necessary to keep small d ^ until the very end, 
even during a calculation of a real part of the polarization 
operator. More formally, the sequences of operations. 



lim 



d^k lim (• 



■ ) and 



lim 

5^0 



d^k lim (• 



(64) 



always lead to a correct static limit, while the sequence 



d^k lim (•••) 

w,(5— 



(65) 



might give in some specific cases an unphysical result 
breaking causality and violating analytic properties of a 
(retarded) response function. 

Let us rigorously prove the statements which have been 
made previously. In studying the limit w of Rell — 
Ren^ + Rell^^ it is sufficient to consider how we approach 
the axis w = from the domains V and A (see Fig. ^ 
that cover the ranges of momenta z G (0, y) and z e 
(y, +cxd), respectively, as w — + 0. In turn, the domains C 
and B in this limit shrink to the points z = and z — y, 
and therefore their consideration is not important for our 
current purpose. 

Let us first consider Rell^ in the domain D at very 
small w. In the range of integration [t2, ti] we can replace 



7^±(T) « 7^l(T) sign(T T wy/z^) - T7^l(T), (66) 
where 



7^l(T) 



22 ^(u;/z)2 -t2' 



(67) 



Meanwhile, in the range of integration [t4, t^] we have 
7^± (r) « 7^2 (r ) = - J- ^z^-{T + y)^. (68) 

2ZT 

We observe that t^+{—) = — r ,_(+) w 1, r+_|_(+) = 

-r_+(-) w w, T++{y) w w/y, T_+(y) w -y; r+_(-) = 
-r__(+) « 1, T__(-) - -r+_(+) « w, r+_(y) « 
-w/y, T — (y) « -y\ r++(0) = -t_+(0) = ^w. Tak- 
ing into account that z < y, we carefully arrange the 
limits of integration, and thus obtain 

- - lim Ren^ = 

z-y r-z-y r-z-y /"^^^^X 

+ / + / +2 / dT7^2(T) 

-y J -y J z-y J -y J 

w I z 

lim I / +1+1 I drTZiir) 




/y J-w/ 



w / z\ 

2 I drTZiir) 

w/y I 



2 / +/ ^^r7^2(T) - lim 4 / drUiir 

\J-y J-y J Jw/y 



2tVz2 - r^dr 2y/f^dr 

zir"^ - y^) 



/y z 



(69) 



by virtue of the table integral ljAll|l . 
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Let us now consider Rell^ in the domain A at very 
small w. Taking into account that now z > y, we obtain 

- - lim Ren^ ^ lim / + / 

mEix{w,w/z) />— min(l,2+y) /'min(l,j2 — y) 



-y -J-y 

ma.x.{w,w/z) p — ma.x.{w ,w / z) />— min(l,2+y) 



2 / -2 / dTlKj) 
"/y -'-y I 



= lim 



rasLyiiw / z) nia'K{w.w/ z) min(l,2H-y) 



'/y 



min(l,2:— y) f- — \/w 



-w/y 



dT2TZ{T). 



(70) 



In the first and the second integrals of the last part of 
(I70|l . we can replace 



7^(r) « 7ti(r) = -1-^^=, (71) 



2z - (u)/z)2' 
while in the third integral 

7^(T) «7^2(T) sign(T). 



(72) 



The fourth and fifth integrals require a more careful con- 
sideration. We note that for small w there holds the 
inequality y/w > w/z, which allows us to neglect w in 



TZ{t). Therefore these two integrals can be rewritten as 



2 hm 



min(l,2:— y) 



y 



c^r7^2(r) 



/>m.in.{l, z—y) 

= 2 1 c^T7^2(T), 

-y 



(73) 



where the integral on the rhs of H73|l is understood in the 
sense of principal value. 
Thus we obtain for z > y 



- lim Ren 

jy w—>0 



max(l,2;) 



2^z2 - y-^dT 



(74) 



z/y Z\/t'^ - 1 

min(l-y..) ^^2 _ ^2^^ fmin{l+y^z) ^^2 _ ^2^^ 







z{T + y) 



z{T~y) 



It is easy to calculate this expression using the table inte- 
gral (|A12(I . The result is expressed in l|A13(l and (|A14p . 

Let us now consider the static limit of Rell-^^ . In Ap- 
pendix^ we show that it equals 



-- lim Ren" = !^V^L!Ze(y - z). (75) 

1/ w^Q z 

We note that this term is nonzero only for z < y (or 
q < 2k a). One can observe that H75|l is exactly canceled 
by the counterterm from Ht)9|) . 

Collecting all contributions and introducing sin^/' = 
y/z for y < z and si'aip± = {1± y)/z for 1 ± y < z, we 
present the static limit of the polarization operator in the 
form 



J 



1 TT 

-Ren(z,0) = 2 9(1 - y - z) + Q{y - \z - 1|) (1 + -sini/;) -2Q{z- 1) arccoshz cosi/' 



Q{z — (1 + /it/)) [ 1 + /iV'M sim/j — cos i/j^ — 2 cos ^ In 



1 -I- zsin(^^ - ^^j) 
2\/2z cos i'f/;^ cos 



(76) 



where Ren(z,0) = limu;_»o Ren(z, w). In what follows 
we also use the abbreviations Ren^(z, 0) and Ren^^(z, 0) 
for lim^j^o Ren^(z, w) and lim„,^o Ren^^(z, w). 



B. Analysis of Eq. iffell 

Let us analyze the expression H76|l . For the values z < 
1 — y we obtain Ren(z,0) = —2^, which ensures the 
fulfillment of the compressibility sum rule. For z > 1 — ?/, 
Ren(2:,0) deviates from the value — 2z/. 

For y = 0, Eq. H7t)|) reproduces the conventional result 



of Stern^ 

-^ReUy=o{z,0) = 2-2^1 - l/z^e{z-l). (77) 

In Fig.2|we compare — Ren(z, 0)/^ for y = 0.1 and y = 0. 
We demonstrate that ~KeIly^o{z,0)/i' is always larger 
than —lielly=o{z,0)/i/, although the difference between 
the two curves (shown in the inset) is quite small. In 
particular, the maximal value of — Ren(z, 0)/i^ at z = 1 
scales with y like 

-iRen(l,0)«2 + ^y3/2^ (78) 
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while at large z the asymptotic behavior of (|76() is 

An important test of our results is provided by the 
Kramers-Kronig relations and the sum rules. For ex- 
ample, using the analytic expressions for Ren(z, 0) and 
Imn(z, w) we have checked the zero-frequency Kramers- 
Kronig relation 



Ren(z, 0) 



dw 



Inin(z, w). 



(79) 



The actual limits of the integration are finite and given by 
the boundaries of the particle-hole continuum. This has 
allowed us to confirm (|79|) numerically with a sufficiently 
good accuracy (the deviation is < 10~^ even for a quite 
simple routine). 

It is instructive to compare Ren(2,0) and Ren'^(2;,0), 



~Ren'^(z,0) = -iRen(z,0) 



-e(y-z). 

(80) 

It is obvious that Rell does not have any anomaly at 
z ~ y (cf. Ref . .12:) ■ while Rell'^ docs have it (see Fig. El . 
One can also observe that Rell^ diverges in the limit 
z ^ and changes the sign at some value of z. Therefore, 
it cannot be regarded itself as a correct static limit. Oth- 
erwise, it would have violated the compressibility sum 
rule and generated an instability of the medium by virtue 
of SO coupling. 

Thus, we have demonstrated that in order to obtain 
the correct static limit of the polarization operator for 
the 2DEG with Rashba SO coupling it is crucially impor- 
tant to follow the thorough definition that implies (|64|l . 
but not H65|l . Otherwise, the contribution Ren^^(z, 0) 
is missing. Being the difference between Ren(z, 0) and 
Ren^(z, 0), it shows up for z < y (or q < 2kji), and dis- 
appears only when kn goes to zero. Therefore, for the 




0.9 0.95 1 1.05 1.1 1.15 z 



FIG. 4: The behavior of -Ren(z, 0)/^ near z = 1 ioi y = 0.1 
(solid line) and y = (dotted line). The difference between 
the two curves is shown in the inset. 
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1 
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;-ReIl^(z,0)/v 







0.05 0.1 0.15 



FIG. 5: The behavior of -Ren{z,0)/iy (solid line) and 
— Ren^(2;,0)/i^ (dashed line) near z = y ior y = 0.1. 



conventional 2DEG (fc/j = 0) there is no difference be- 
tween Hq Q and Ilq^o [and between H64I) and (|65|l as well]. 
So the property H63() is a peculiar feature of the 2DEG 
with kfi 0. We suppose that it might be related to the 
singularity of the spectrum ((2Jl at fc = 0. 



C. Effective interaction: modification of the Friedel 
2fci7-oscillations 

Let us conveniently rewrite the RPA expression ^ for 
the dielectric function in terms of f. 



2V2' 



e{z, w) — 1 -n(z, w), 



(81) 



where 



_ \/2n 



is the 2D Wigncr-Scitz parameter. 



The latter controls the accuracy of RPA, which becomes 
better with decreasing r^. In the limit w ^ 0, Eq. H81|) 
describes the static screening of the Coulomb interaction. 
It is well-known^ that the singular behavior of the deriva- 
tive. 



1 d 

V dz 



ny=o(^,o)|^^^_^^ 



(82) 



a,t z = 1 + a with small a > gives rise to the Friedel 
oscillations of a screening potential, their leading asymp- 
totic term being 



2e^kF 2fse 



-2kpd 



sin 2kFP 



[£,.0(1,0)]2 {2kFpf 



(83) 



where p = \/ x\ + and d is a distance from a probe 
charge to the 2DEG plane. The SO modification to 
follows from 



V dz 



n(z,o) 



(84) 
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and results in 



Voscip) 



where the factor 



Q{y) = 



=o(l,0) 



1 2 



e(l,0) 



(85) 



(86) 



is always smaller than 1. It means that due to the SO 
coupling the amplitude of the Friedel oscillations is di- 
minished (cf . Ref . 112} , although the amount of such de- 
crease is quite small (~ 0.5% for y — 0.07 and = 0.2). 
We note that due to (|78|l it is possible to approximate 



1 2 



£(1,0) 



1 



3(1 + 2f,)' 



(87) 



Rigourously speaking, it is necessary to take into account 
in l|H5|l the dependence oi kp on as well. 

The second derivative j^Il{z, 0) also diverges at the 
points z = liby-l-0"'", thus contributing to the sublcading 
asymptotic terms of the oscillating potential^—. 



V. STRUCTURE FACTOR AND SO-INDUCED 
DAMPING OF PLASMONS 

An important quantity which cither can be directly 
observed or enters into expressions for other observable 
quantities is the structure factor^ 



S{z, w) = — Im [l/e{z, w) 



It depends on the Wigner-Seitz parameter and con- 
tains information about both particle-hole excitations 
and collective excitations (plasmons). The spectrum of 
the latter is found from the equation Ree — 0. For 
the conventional 2DEG it can be derived on the basis 
of Ref. "Si and equals^®- 



z{z + 2Ps) /4f2+4f5z3 + z4 ^ 
wpliz) = ^— '-J^——— e(z* - z), (89) 



2r, 



z{z + 4fs) 



where the endpoint of the spectrum z* is the real positive 
root of the equation z'^{z+Ars) = 4f^. Provided Ime — 0, 
the plasmon spectrum is undamped, and the structure 
factor is (S-peaked at w = Wpi{z), 



Spi{z,w) = a{z)6{w - Wpi{z)), 



(90) 



with the weight factor a(z) — 7r[Re9e/9?i;|^^^,^j(2)] 
given by 



a{z) 



7rVI[16ff-z4(z-f-4f,)^] 
2f5^(4f2 -I- 4fsz3 + z^)(z + 4rs)3 ' 



(91) 



The plasmon spectrum can be visualized on a contour 
plot of S{z, w) by adding an artificial infinitesimal damp- 
ing S to Ime. It becomes very helpful when an exact 
analytic expression similar to (|89|) is not known. 




FIG. 6: Contour plots of S{z,w) showing the SO-induced 
wedge-shaped damping region (bounded by the dashed lines). 
The plasmon mode is depicted by the bold line. Insets show 
the cross-sections S{z = 0.1, w) as a function of w. Left panel: 
y = 0.07, rs = 0.2. Right panel: i/ = 0.04, r, = 0.6. 



In the presence of SO coupling, the structure factor 
depends on y as well. Solving approximately the equation 
Re£ = at 2/ 7^ 0, we can establish how the plasmon 
spectrum is modified by SO coupling. In fact, we observe 
that the SO modification of (|89() is quite small. Later on, 
we will comment on how small it is, and explain why this 
effect is not really important (cf. Ref. 15). 

A much more important effect is a damping of plas- 
mons generated by SO coupling. As it has been already 
discussed, SO coupling extends the continuum of the 
particle-hole excitations up to the new boundaries shown 
in Fig. 12 Therefore, the plasmon spectrum is expected 
to acquire a finite width, whenever it enters into the SO- 
induced region of damping. 

In Fig.Elwe show the contour-plots of the structure fac- 
tor S{z, w) depicting the plasmon spectrum by the bold 
line, in black, where Ime = (undamped plasmon, the 
structure factor is J-peaked) and in gray, where Ime ^ 
(SO-damped plasmon, the structure factor has a finite 
height and width). Depending on the values of Vg and y, 
the two different cases are possible: (I) the plasmon en- 
ters only once into the SO-induced damping region (left 
panel); (II) it enters twice, escaping for a while after the 
first entrance (right panel). 

Within the conventional boundaries w = z'^ zL z of the 
particle-hole continuum, the structure factor S{z, w) is 
modified by SO coupling only slightly and can be ap- 
proximated by the conventional expressions^. 

In the SO-induced region of damping, S{z,w) is very 
well approximated by the Lorentzian function describing 
the SO-damped plasmon with the width 7(2), 



S{z, w) so-damp pi 



a{z) 



7(z) 



TT {W - Wpl{z))'^ +j'^{z) 



, (92) 



where Wpilz) and a{z) are supposed to be practically 
independent of y and therefore given by l|5^ and 
while 



7(2;) = I!iflinre(z, Wp;(2)) 



(93) 
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FIG. 7; The cross section S{z — 0.1, w) (solid line) compared 
to the approximation S{z — 0.1,11;) so -damp pi (dashed line). 
Parameters: i/ = 0.07, = 0.2. 



essentially depends on y via Ime ^ 0. In Fig. [7| we 
present the enlarged plot with the cross section S{z — 
0.1, ■«;) from the inset of the left panel of Fig. Inland com- 
pare it with the approximate S {0.1, w)so-damp pi given 
by (|92|l . The inset of Fig. [7| shows both curves on a more 
fine scale. Comparing the positions of their peaks, we 
conclude that the shift of the plasmon dispersion due to 
SO coupling from Wpi (z) (|89(l is one or two orders of mag- 
nitude smaller than 7(2:) [unless 7(2;) — 0], and therefore 
can hardly be resolved experimentally. The almost equal 
height of the peaks confirms that (|91|l is also a very good 
approximation for the weight factor at y 7^ 0. 

Let us find the values of y and Ts that correspond 
to the typical cases (I) and (II) shown in the left and 
right panels of Fig. (BJ respectively. For this purpose, it 
is necessary to establish when the curves wi{z) (|49|l and 
Wpi(z) 189(1 touch each other. The result is represented 
in Fig. |S1 where the plane (y,rs) is divided into the cor- 
responding domains I and II. Changing y and Vg , we can 
tune the relative position of the SO-induced particle-hole 
continuum and the plasmon spectrum. The insets show 
7(2) = 7(2) X 10* for the values of parameters the same as 
in Fig.El (I) y = 0.07, rs = 0.2; (II) y = 0.04, r, = 0.6. 
The function 7(2) is defined for z < z* and vanishes 
where the plasmon is undamped. 

A direct measurement of SO-induced plasmon width 
can be provided by means of inelastic light (Raman) 
scattering (see, e.g., Ref. [T^ . Experimentally it is 
possible to measure the structure factor in the range 
q = — 2 X 10^m~^ which corresponds to the range of 
z = — 0.1 (for rg ~ 1). Keeping z at a fixed value 
and varying y and r^, one should observe different pat- 
terns of S{z,w) with either damped or undamped plas- 
mon, like those shown in the insets to Fig. El where the 
sections of S{z,w) at constant z = 0.1 are presented. 
In order to estimate realistic parameters, we focus on 
InAs-based 2DEG, where the strength of Rashba cou- 
pling an has quite large values up to ~ 3 x 10~^^eVm 
(Ref. 59) and prevails over the Dresselhaus term2i. Most 
of experiments deal with the 2DEG's densities ranging 
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FIG. 8: The plane {y, r^) is divided into the domains I and II, 
where the plasmon has one or two undamped pieces, respec- 
tively (cf. Fig. EJ. Insets show SO-induced plasmon width 
7(2) = 7(2:) X 10* for the same parameters as in Fig. |S| 



from n = 0.7 X lO^^m--^ (Ref. |2a) to n 
(Ref. [13). Taking ~ 12 and m* w O.OSrrie, we obtain 
the range of y from 0.04 to 0.075 and the range of 
from 0.38 to 0.3. These values belong to the domain I 
in Fig. |S1 where the SO-damping of the plasmon is es- 
pecially important. Choosing the value of 7 ~ 10~*, 
we establish that 2Tpi = h/{A'yEp) is of the order of 10 
ps. This can, in principle, be resolved experimentally, 
since the Raman measurements done in II- VI quantum 
wells^ give a finite lifetime of the plasmon mode with 
the typical value 2rp/ ~ 0.3 ps. We also point out that 
the II- VI compounds (e.g., HgTe) represent an even more 
promising candidate for observation of the discussed ef- 
fects since the reported SO coupling strength is compar- 
atively large^. 

Plasmon broadening may be also caused by ther- 
mal effects. We use the results for 2DEG at finite 
temperature^^ in order to estimate the relative contribu- 
tions of SO-induced and thermal damping. The former 
dominates over the latter below some characteristic tem- 
perature that is about 90 K for the parameters of Ref. I20I 
This estimate is made at constant z — 0.1. Note that the 
thermal damping strongly depends on the values of 
and temperature, whereas the SO-induced damping has 
quite a weak dependence (e.g., compare the insets in 
Fig. |SJ). This may provide a guide for an experimental 
separation of these two sources of broadening. 

A finite lifetime of plasmons would modify the energy- 
loss function. The energy loss per unit length for a par- 
ticle moving toward a plane with 2DEG with velocity v 
is given by 



^ = Cm Hzdz 
ax \ V / Jq 



'"^'"^ uS{z,uz)du 



(94) 



where C = l&ir'^e^n/ Sao- In the presence of SO coupling, 
the most essential modification is expected in the plas- 
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FIG. 9: The plasmon contribution to the energy-loss function 
at ?/ = (dashed hne) compared with the SO-damped plas- 
mon contribution at y — 0.07 (sohd hne). For the both curves 
Ts = 0.2. 

moil sector. In Fig.l^we plot the plasmon contributions 
{—dW/dx)pi for y = and y = 0.07 as a function of 
v/vp. Both curves are calculated at = 0.2 and nor- 
malized by the peak value of {—dW/dx)pi at y = 0. We 
observe that the sharp peak is smoothed down and the 
steplike behavior becomes more smeared with increasing 
SO couphng strength. 
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FIG. 10: SO-induced damping of the LO phonon 7^0(2) = 
7^0(2) X 10"' for lOo = 0.07, = 0.4, Soo = 12, and eo = 
15. Solid and dashed lines correspond to the cases y — 0.07 
(jwo — y\ < y^) and y — 0.06 (|wo — y\ > y^), respectively. 
Insets show the location of the phonon spectrum (solid line) 
relative to the SO-induced damping region (dashed lines) and 
to the strip y ~ y^ < w < y + y^ (dotted lines) for y = 0.07 
(left) and y = 0.06 (right). 



VI. SO-INDUCED DAMPING OF LO PHONONS 

Our results for the structure factor allow us to predict 
that the other collective mode, longitudinal optical (LO) 
phonon, will also experience a pronounced SO-induced 
damping. Its dispersion and lifetime can be obtained 
from the renormalized propagator that is expressed in 
RPA as26 



Diq,Lu) 



"iuJLO 



1 



(95) 



where — ^KVqUiLOj ^^nd k — 1 — £00 1 £0 is a material 
parameter. On the basis of this expression, we can estab- 
lish that the LO-phonon width (normalized by 2kp/'m*) 
equals 



KS{z,w){wo'Re e) 



2w(Ree)'^ + KWQRe{de / dw) 



(96) 



where wq = (m* /2kp) ujlo and Wph{z) is the renormal- 
ized phonon spectrum. Like in the case of plasmons, 
we may neglect the dependence of Wph{z) on SO cou- 
pling, and find the phonon spectrum from the equation 
w'^ — Wg[l + ^((Ree)^-'^ — 1)], where Ree is taken at 
2/ = 0. 

The lifetime effects of the LO phonons can be very im- 
portant for a coupled dynamics of carriers and phonons. 
An interesting possibility arises when wlo ^aakp (or 



Wq ^ y). This, for example, holds in InAs, where the 
value ujlo ~ 28 meV gives wq « 0.07. Changing the 
Rashba coupling strength by an applied electric field, 
one can manipulate the lifetime of an optical phonon, 
which would result in a modification of transport prop- 
erties by virtue of the electron-phonon coupling. 



In Fig.^lwe show the function 7^0(2) = Jlo{z) x 10^. 
In our calculations we used the following parameters: 
Wo = 0.07, Ts = 0.4, £00 = 12, and £0 = 15 (for InAs). 
Choosing the value jlo = 10"'', we find that for LO 
phonons the SO-induced lifetime 2tlo is of the order of 
10 ps as well. We note that the typical lifetime for the 
LO phonons measured in AlAs and GaAs by the Raman 
spectroscopy is of the same order—. One can observe that 
for y = 0.07 {\wo — y\ < y^) the damping is always finite 
(solid line), while for y — 0.06 {\wq — y\> y"^) it becomes 
nonzero only after some value z w 0.007 (dashed line). In 
the two insets we show the relative location of the phonon 
spectrum (solid line) and of the SO-induced damping re- 
gion (dashed lines). We also depict the boundaries of the 
strip w — y zL y'^ (dotted lines), and one can see that 
for y — 0.07 the phonon spectrum leaves it at the value 
z « 0.017 (left inset), while for y = 0.06 the phonon 
spectrum never gets inside the strip (right inset). This 
indicates that the the approximation H59|l is insufficient 
for describing the discussed effect. 
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VII. CONCLUSION 

We have calculated the dielectric function of 2DEG 
with Rashba SO coupling in RPA at zero temperature. 
We have described the new features of screening that 
appear due to SO coupling. In particular, we have dis- 
cussed in detail the extension of the region of the particle- 
hole excitations, which leads, in turn, to the additional 
broadening of the plasmon mode. The same mechanism 
generates the SO-induced lifetime of the longitudinal op- 
tical phonons. Speaking generally, we have seen how SO 
coupling tends to suppress collective excitations due to 
the relaxation via intersubband transitions. At the same 
time, SO coupling does not affect much the position of 
the collective mode dispersions, and therefore this effect 
is less important. 

Another conclusion that can be derived from our stud- 
ies is that the 2DEG with Rashba SO coupling requires 
a more careful treatment as far as various approxima- 
tions are concerned. Usually they fall into the follow- 
ing categories: (1) the so-called ^-approximation based 
on the linearization of the spectrum; (2) an expansion 
of physical quantities in powers of SO-coupling strength 
before their evaluation; (3) reshuffling the operations of 
the momentum integration and taking the limit of zero 
frequency. The examples that indicate that none of these 
approximations provide a fully reliable result for the sys- 
tem in question are the following: (1) being a version of 
the ^-approximation, Eq. (|56|l has a very limited range 



of applicability; (2) the static polarization operator at 
small y cannot be obtained by means of the power ex- 
pansion in y, because its maximal value scales as y^^"^ 
[see Eq. (I78f) ]: (3) reshuffling the order of operations leads 
to the unphysical result with violated analytic properties. 
Although the underlying reason of the failure is not quite 
clear, the above examples warn against the blind usage 
of these popular approximations in this system. 

Finally, the comparison of our theoretical estimates for 
SO-induced lifetimes of plasmon and LO-phonons with 
the values of the recent Raman scattering measurements 
suggests that the effects described above can be observed 
exp er imcnt ally. 
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APPENDIX A: SOME TECHNICAL DETAILS OF 
THE DERIVATION OF n(q, uj) 

Intermediate steps that lead from from HlOf) to Ullf) - 
(|14|l are to express 



A 



g,{v,z,w,y) = 7 / d(j3 



2tt 



1 



V + 2zx 



1 - 



|v + 2z 

V + 2zx 



S{2Xw — 2z'^ — 2vzx + fivy — /J.2/|v -I- 2z| 



5{2Xw - 2z - 2vzx + fivy + ^y|v + 2z|) 



fi{v,z,w,y) 



|v + 2z|, 

vzx + z"^ — ^'Uizx + 7j) — \w 



TT {2vzx + 2z'^ — fiyv — 2Xw)^ — y'^(v'^ + Az^ + Avzx) 



(Al) 
(A2) 



The change of variables H34() implies 

t'^ + yr — Xw 
V = , 



XL -I = - 



t'^ + yr ~ Xw ' 
(r^ + yr — Xw)^ ' 



(A3) 



'^dv dr 

= — (A,+ - A,_), (A4) 

V yr 

TZ-Xw{T + y)/z 



(A6) 



-|- yr — Xw 

where Xia- are defined in (|16|l . and Ci and t^. are given by 
((T^ and linZI), respectively. 



In this Appendix we also quote table integrals that 
we use in our calculations. They can be either found in 
Ref.|l3lor elaborated on its basis. We quote the standard 
integral (jA8|) containing the (5-function as well. 



d(j}5{a + cos ( 



=6(1 -H), 



27r 



a + cos(t> Va^ - 1 



(A8) 



e(|a| - l)sign(a), (A9) 
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X2 + /Jx + 7 ^(l+7)2-/?2 



(AlO) 



(J+l)Vl-/3 + 7+(^-l)Vl + /3+J7 (4^>^2) 



dr Va^ - 1 
T - a Vl - r2 



(All) 



2 arctan - 



-, (|r|<l<|a|) 



dr VI -a2 
T - a Vl - r2 



In 



(a — T + aVl — — T 



2a(a-T)(l + Vl-T2) 
Using l)A12|l it is easy to establish that 

.max(l.^) 2^^2 _ y2^^ 
Jz/y Z\/t2 - 1 



(A12) 
(|rUa|<l). 



(A13) 



9(z — l)arccoshz — arccosh— 

y. 



and 



z(r±?/) 



= i-e(iTy-^) (A14) 



1 - ^ In - 



+ \Jz^ -y , Try 



± 



2z 



e(z-(iTz/)) 



1 — ^ ± — arcsm 



2^/1- ^In- 



(It;;)'' 



2y l + \/i- 



APPENDIX B: AN ALTERNATIVE 
REPRESENTATION OF fo,]^ 

For > jif {i = 3, 4) the roots of the equation x"^ + 
(3iX + 7i = are complex. Omitting the index i (in order 
to avoid confusion with imaginary i), we write them in 
the form [a = ±) 



A. = --+^ay7-_, 



(Bl) 



and note that A_ = A!j_, |Acr|2 = 7. The functions /g 4 
/ are then represented by 



/ 



iC v-^ . , . \a + 5 



dz 



dz 



'\z\ = l Z 

where the secondary roots {a' = ±) 



(B2) 



W = A, + (t'VA2~T (B3) 

obey the equation z-^^, — 2\crZcra' + 1 = 0. Note the 
properties z-cr' = z^^^i and Za+Za-^ — 1. We also denote 



where a and b are real. 



Checking which of the roots Zo-o-' lie inside the unit 
circle \z\ — 1, we find an expression for IIB2|l . 



/ = - 



C sign(a) 



A, 



A_ 



Ai -1 



XI - 1 



C sign(a) _^ X+ + S 
-Im- 



ImAa 



(B4) 



One can notice that it resembles its /■'^-counterpart H21|l , 
and therefore a (complex) change of variables similar to 
(|34|l would lead to the form of integrand the same as in 
(|36|l . However, the corresponding mapping of the inte- 
gration contour would assume that an integration vari- 
able T runs in the complex plane along the circle of the 
radius \/w. This would probably require an analytic 
continuation of the elliptic functions (see Appendix . 
which makes very cumbersome an explicit representation 
of 



APPENDIX C: INTEGRALS EXPRESSED IN 
TERMS OF ELLIPTIC FUNCTIONS 



Let us consider the integral 



. {x- Xa){x - Xb) 

1=1 — dx, 



(CI) 



assuming that X4 < x^ < X2 < xi. We introduce the 
notations Xrs = x^ — x^ and 

= (C2) 

X'ilXi2 

Let us consider the intervals where the polynomial 
Y\f^^i{x—Xk) is positively defined. For a; < ^4 and xi < x 



the integral IjCip equals 



1 



' {x — Xi){x — Xz){x — X4) 
X ~ X2 



(C3) 



^/X3lXi2 



[(2Xa2Xb2 + X42X2l)F{(p, k) - X'iiXi2E{Lp, k) 



Xil 



2^21 (2Xa2 + 2X62 " 2:32 " X42 -|-X2l)n(^, ,/c)], 

2:42 
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where 



ip = arcsm 



XA2 X — X\ 
X41 X — X2' 



(C4) 



and the sign "+" has to be chosen for xi < x, while the 
sign "— " has to be chosen for x < X4. 
For X3 < X < X2 



I = 



{Xi - X){X2 - x){x - X3) 



X ~ X4, 



(C5) 



1 



[{2XaiXb4 - X42X4:3)F{(p, k) - X3iX42E{(p, k) 



X32 



+ X43(2Xa4 + 2XM + 2^43 + Xi2 + a;4l)n(<^, , /c)], 

X42 

where 



(f = arcsm 



0:42 X- X3 
X32 X — X4 



The eUiptic functions 



F{cp,k) = 



da 



y/l- k'^si 

V 



2 

sm a 



E{(p,k) = I da\/ 1 - fc2 sin^ a , 
Jo 

U{Lp,n,k) = / „ " , == 

"'o (1 — nsin ajyl — k-'sm a 



(C6) 

(C7) 
(C8) 
,(C9) 



are defined according to Ref. 



APPENDIX D: STATIC LIMIT OF ReH" 

We start out from the expressions (|22|) and H29I) . For 
smaU w < 1^, we need to consider only the contribu- 
tion from the function /{^ specified by /i = — , i.e., 

lim Ren^^ = lim / vfi'{v,z,w,y)dv. (Dl) 

Let us introduce a new integration variable f ~ tt-^ 

^ ^ 2^w 

m ||DT|| . Thus, we make the limits of integration inde- 
pendent of w, and at this stage it becomes possible to 



exchange the operation lim^^o a-nd the integration over 



- - lim Ren 

V w—>0 



II 



lim / d^4>{C,z,w,y) 



0o(C,2;,?/)d^, 



(D2) 



where 



X fi'{y + 2^V^,z,w,y), (D3) 
(f>t){Lz,y) = lim (/)(^,z,w,?/). (D4) 
We observe that 



^^^\z-y\ 

z'^ + {z - y f - ^^y'^ + 4z{z - y)) 
2z2(z-y)2 

^/Q^=V^Az + y) 
7? + {z + yf - (^{y^ + Az{z + y)) 



1 - 



z{z-y) 



,(D5) 



2z^{z + yf 



z{z + y) 



For z > y, 



2 



(D6) 



(D7) 



(/f^+v/Q^) -Vz « Az{z^~y^), (D8) 
and therefore 4>{){(,, z, y) ~ 0. For z < y, 



and therefore 



z{y^ — z^) 



0o(C,^,y) 



When inserted in HD2|) . it yields (|75|l . 



(D9) 
u;, (DIO) 
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